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ABSTRACT 


The  problem  of  obtaining  a  minimal  Gaussian  basis 
set  is  treated  in  the  present  thesis.  There  are  a  number 
of  such  sets  already  available,  but  they  all  have  serious 
limitations  and  shortcomings.  The  most  celebrated  of 
these  existing  minimal  Gaussian  basis  sets  are  reviewed, 
and  their  defects  are  found  and  proved.  A  review  of  the 
method  of  obtaining  basis  set  expansions  is  also  made. 

A  new  basis  set  is  proposed  here.  Atomic  orbitals 
are  expanded  in  terms  of  three  Gaussian  functions.  The 
expansions  are  determined  by  a  least  squares  fitting 
process  developed  by  S.  Huzinaga  et  al.  The  orbital 
exponents  and  expansion  coefficients  are  then  energetically 
optimized . 

Applications  of  the  new  basis  set  are  made  on  a 
number  of  atoms.  The  results  are  very  encouraging 
considering  the  level  of  sophistication  of  the  method. 

It  is  demonstrated  that  STO-SZ  and  ST0-3G  are  very  in¬ 
effective  when  applied  to  third  row  atoms.  In  contrast, 
the  values  produced  by  this  work  for  the  mentioned  atoms 
are  very  good.  The  method,  therefore,  proves  highly 
appropriate  for  treating  transition  metal  elements. 

Quantum  mechanical  calculations  of  some  closed-shell 
molecular  systems  were  also  performed.  The  molecular 
and  molecular  orbital  energy  values  thus  determined  were 
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found  to  be  significantly  better  than  those  obtained 
using  ST0-3G.  Molecular  geometries  are  also  calculated 
to  a  satisfactory  degree  of  accuracy. 

The  method  leads  to  contracted  Gaussian  basis  sets, 
which  are  compact,  easy  to  use,  and  very  economical  in 
terms  of  computing  time.  It  can  thus  safely  be  recom¬ 
mended  for  calculations  of  polyatomic  molecules.  It 
should  particularly  be  used  in  cases  involving  third  row 
or  heavier  atoms,  where  the  STO-KG  method  breaks  down. 
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I .  INTRODUCTION 

The  first  steps  in  the  history  of  quantum  mechanical 
calculations  on  the  electronic  structure  of  molecules 
aimed  at  the  determination  of  a  wave  function  for  a  single 
given  arrangement  of  the  atoms  in  space.  In  most  cases 
this  arrangement  was  the  equilibrium  geometry,  as 
determined  experimentally. 

The  development  of  high  speed  electronic  computers 
has  been  extremely  rapid  during  the  last  few  decades. 

The  resulting  emphasis  on  quantum  mechanical  calculations 
led  to  the  refinement  of  the  theories  and  methods  used 
in  such  calculations.  Potential  energy  surfaces  with 
respect  to  geometrical  parameters  have  consequently  been 
extensively  explored,  in  order  to  locate  the  minima  of 
these  surfaces,  and  thus  establish  the  equilibrium 
values  of  the  structural  parameters. 

The  search  for  potential  energy  minima  can  be  done 
in  two  ways.  One  can  assume  that  the  values  of  a  number 
of  parameters  are  constant  and  vary  the  rest  until  a 
minimum  is  achieved.  Alternatively,  the  energy  can  be 
minimized  with  a  complete  and  simultaneous  variation  of 
all  parameters.  It  is  obvious  that  the  latter  method 
would  be  preferable,  whenever  such  a  process  is  possible. 

One  is  thus  led  to  an  a  priori  prediction  of 
molecular  geometries.  In  the  process,  the  only  experimental 
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data  used  are  those  relating  to  the  values  of  the  fund¬ 
amental  physical  or  chemical  constants.  This  means  that 
the  calculated  geometry  is  purely  a  result  of  the  theory 
and  method.  The  theoretical  structures  can  then  be  used 
in  two  ways : 

(a)  to  provide  a  systematic  evaluation  of  the  capabilities 
and  potentials  of  the  theory  and  method.  Also,  to 
uncover  any  possible  defects  or  inherent  limitations 
of  the  theory,  at  its  certain  level  of  sophistication. 

(b)  Once  the  theory  and  method  have  been  successfully 
tested,  they  can  be  used  to  calculate  the  structures 
of  molecules  that  have  not  so  far  been  determined  by 
experiment.  These  calculations  can  only  be  good  to 
the  established  level  of  accuracy  of  the  theory  and 
method. 

Most  of  the  theoretical  structural  predictions  to 
date  have  been  carried  out  at  the  Hartree-Fock  level, 
using  a  single  determinant  wave  function,  with  a  basis 
set  expansion.  A  basis  set  expansion  is  usually  a  linear 
combination  of  functions  that  represents  the  behaviour  of 
an  atomic  or  molecular  orbital.  The  accuracy  of  the 
results  is  closely  related  to  the  basis  set  expansion. 

For  a  large  basis  set,  the  results  will  be  independent  of 
detail  and  characteristic  of  a  geometry  at  the  Hartree- 
Fock  limit.  For  a  small  basis  set,  the  results  strongly 
depend  on  the  basis  set  employed. 
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Since  in  the  majority  of  cases  we  are  forced  to  use 
a  limited  basis  set,  the  importance  of  its  choice  is 
obviously  critical  to  the  calculations;  hence  it 
should  be  made  with  the  utmost  care.  Physical  and 
chemical  intuition  and  insight  are,  therefore,  very 
important . 

In  the  next  chapter  we  shall  demonstrate  the  use  of 
two  different  kinds  of  basis  functions  in  simple  quantum 


mechanical  calculations. 
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II  STO's  AND  GTO's 


II.  1  The  Hydrogen  Molecule  with  Slater-Type  Functions. 

Because  of  its  simplicity,  the  hydrogen  molecule  has 
been  heavily  used  as  a  prototype  in  demonstrations  of 
various  theories  and  methods.  Although  some  of  its 
properties  are  quite  unique,  it  can  still  provide  a  useful 
example  from  which  ideas  about  the  use  of  single¬ 
electron  wave  functions  in  the  expansion  of  molecular 
wave  functions  may  be  drawn. 

Consider  the  hydrogen  molecule  at  infinite  inter- 
nuclear  separation.  In  this  case,  the  two  atoms  can  be 
considered  independent  of  each  other.  Then,  the  wave 
function  of  the  ground  state  of  the  molecule  is'*': 


-6(l)a(2)| 


(1) 


where 


and  N  is  the  normalization  factor. 

From  equations  (1)  and  (2) ,  it  follows  that  functions 
of  the  form: 
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rn  le  rjrYlm(0,<)>)  <3> 

are  suitable  for  expanding  the  molecular  wave  function. 

.  2-5 

These  functions  (3) ,  are  called  Slater-type  functions 
It  must  be  remembered,  however,  that  the  above 
discussion  only  applies  at  infinite  internuclear  distance. 
When  the  atoms  approach  each  other  for  molecular  formation, 
the  electronic  distribution  is  bound  to  change.  This 
means  that  the  shape  of  the  orbitals  will  be  altered. 

For  finite  internuclear  distances,  therefore,  the  wave 
function  (1)  is  still  a  good  approximation  to  the  true 
molecular  wave  function,  if  it  is  adjusted  to  accommodate 
for  the  orbital  contraction  that  takes  place  upon  molecular 
formation . 

The  experimental  dissociation  energy  for  the  hydrogen 
molecule  is: 

DexP  =  4.72  eV 
e 

The  dissociation  energy  obtained  using  function  (1)  is: 

D  =  3.14  eV  with  C  =  1.0  in  (2) 
e 

To  account  for  the  change  of  atomic  orbitals,  upon 
molecular  formation,  minimize  the  energy  with  respect  to 
the  orbital  exponent  £,  in  formula  (3) .  This  leads  to 
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the  following  values: 

C  =  1.17,  De  =  3.76  eV 

The  same  result  could  be  attained  by  a  different 
method.  Instead  of  using  one  Slater-type  function,  use 
the  linear  combination  of  two  such  functions  on  each 
centre,  i.e. 


(4) 


ls±  =  Cilsi^i^  +  C2lsi^2^ 


There  is  an  advantage  in  this  latter  method,  because  one 
substitutes  the  optimization  of  the  non-linear  parameter 
C  in  (3) ,  by  the  optimization  of  two  linear  parameters 
(C  , )  in  (4),  which  is  an  easier  task.  Choosing  =  1.0 
and  =  1.5,  one  obtains: 


D  =  3.74  eV 
e 


C  =  0.663,  C2  =  0.351 


It  is  possible  to  further  improve  the  wave  function. 
This  is  achieved  by  taking  into  account  the  polarization 
effects,  with  the  inclusion  of  the  a(2p)  function  to 
provide  a  better  buildup  of  charge  in  the  bonding  region. 
The  dissociation  energy  thus  attained  is^: 
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D  =  4.02  eV 
e 

It  is  also  possible  to  consider  a  completely  general 
function  '  ,  and  determine  its  most  appropriate  form. 

The  Heitler-London  valence-bond  method  we  just 
described  depended  on  the  assumption  that  a  wave  function , 
valid  in  the  separated  atom  region,  could  be  extended  to 
apply  to  the  molecular  region.  But  it  is  equally  possible 
to  start  at  the  united-atom  region  and  extend  the  wave 
function  into  the  molecular  region.  This  leads  to  the 
molecular-orbital  method,  which  is  the  complete  analogue 
of  the  atomic-orbital  method  for  atoms.  It  must  be  noted, 
however,  that  molecular  orbitals  are  normally  polycentric, 
as  compared  to  monocentric  atomic  orbitals. 

Let  us  consider  the  molecular  orbital  for  the  ground 
state  of  the  hydrogen  molecule.  It  is  reasonable  to 
assume  that  the  wave  function  near  nucleus  A  will  resemble 
the  atomic  orbital  lsA  of  that  atom.  The  same  argument 
applies  for  the  wave  function  in  the  vicinity  of  nucleus 
B.  The  Ritz  principle  thus  leads  us  to  write: 

»+  =  C1  1SA  +  °2  1SB  Ua) 

where  the  Is  function  is  the  same  as  in  (2).  The  co¬ 
efficients  C  and  C2  are  variable  parameters.  In  this 
,  due  to  the  symmetry  of  the  hydrogen  molecule  with 
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respect  to  interchange  between  A  and  B,  we  can  write: 

¥  ,  =  Is,  +  lsn  (lb) 

*r  A  d 

Evaluating  the  dissociation  energy  as  before,  with  the 
value  of  £  in  (2)  equal  to  unity,  we  obtain: 

D  =  2.70  eV 
e 

However,  this  value  can  be  considerably  improved  if  we 
optimize  the  orbital  exponent  C.  In  this  case: 

D  =  3.14  eV 
e 

Further  refinements  to  these  values  can  also  be  made,  as 
they  were  made  on  the  Heitler-London  wave  function. 

In  the  next  section  the  use  of  a  different  kind  of 
function  will  be  demonstrated  through  its  application  to 
the  hydrogen  atom. 

II. 2  The  Hydrogen  Atom  with  Gaussian-Type  Functions. 

In  1950,  S.F.  Boys  first  demonstrated  the  usefulness 
of  Gaussian-type  functions  in  quantum  mechanical  cal¬ 
culations.  The  general  form  of  these  functions  is: 

i  2 

1  m  n  -ar 
<f>  =  x  y  z  e 


(5) 
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Let  us  employ  one  such  function  to  describe  the 
ground  state  of  the  hydrogen  atom.  Optimizing  the 
Gaussian-orbital  exponent  in  (5),  we  obtain: 

a  =  0.283  E  =  *-0.42441  a.u. 

The  atomic  energy  of  hydrogen  is  known  to  be 
-0.5  a.u.  Thus ,  use  of  only  one  Gaussian  function  yields 
approximately  85%  of  the  energy.  If,  however,  we  use  a 
linear  combination  of  four  such  functions,  i.e. 


4 

=  I 


i=l 


C.cj>. 

iYi 


(6) 


and  optimize  both  the  orbital  exponents  and  the  expansion 
coefficients,  we  have: 

E  =  -0.49928  a.u.10 

We  thus  retrieve  99.86%  of  the  total  atomic  energy. 

For  reasons  of  simplicity,  only  the  hydrogen  atom 
was  discussed  above  using  Gaussian  functions.  A  treatment 
of  the  molecular  case  would  follow  along  the  lines  of  the 
previous  section,  the  only  difference  being  that 
Gaussian-type  functions  would  now  be  used. 
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There  is  one  unavoidable  conclusion  from  the  preceding 
example.  It  is  possible  to  obtain  an  accurate  enough 
description  of  an  atom  provided  that  a  sufficient  number 
of  Gaussian-type  functions  is  used. 

In  the  following  section  we  attempt  to  make  a 
comparison  between  the  two  kinds  of  functions  we  have 
considered  thus  far.  Some  of  their  properties  will  be 
mentioned  and  reasons  for  the  introduction  of  Gaussian 
functions  to  molecular  calculations  will  be  given. 

II. 3  Slater-Type  Functions  versus  Gaussian-Type  Functions. 

It  has  already  been  demonstrated  that,  at  least  in 
principle,  Slater-type  functions  are  highly  appropriate 
for  use  in  molecular  calculations.  They  have  indeed  been 
extensively  used  in  such  calculations  with  significant 
success.  It  is  essential  to  mention  that  due  to  their 
exponential  character  they  provide  a  very  accurate 
description  of  the  behaviour  of  the  electrons  near  the 
nucleus.  This  is  a  significant  property  that  Gaussian- 
type  functions  do  not  possess.  It  is  unfortunate, 
however,  that  Slater  functions  lead  to  extremely  difficult 
many— centre  two— electron  integrals,  the  solution  of  which, 
when  possible,  requires  unreasonably  long  computer  time. 
The  only  exception  is  the  case  of  diatomic  molecules,  which 
have  been  calculated  using  Slater-type  functions.  The 
results  of  these  calculations  have  been  extremely 
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successful,  and  they  are  the  best  theoretical  calculations 
for  diatomic  molecules  available. 

9 

On  the  other  hand,  the  work  of  S.F.  Boys  has  proven 
that  molecular  integrals  can  be  evaluated  with  relative 
simplicity,  if  Gaussian  functions  are  used.  But  we  have 
already  seen  that  Gaussians  are  inherently  poorer  than 
STO's  in  describing  atomic  and  molecular  systems.  It 
is  necessary,  therefore,  to  use  a  larger  number  of  the 
former  for  results  of  comparable  accuracy. 

In  Table  1,  the  convergence  of  the  Hartree-Fock 
energy  of  the  oxygen  atom  as  a  function  of  the  number  of 
Gaussian  functions  used  in  the  calculation  is  shown"^. 

An  examination  of  the  table  verifies  two  points  that  have 
been  made  above : 

(a)  that  the  energy  convergence  improves  drastically  with 
the  increasing  number  of  Gaussians  used. 

(b)  that  Gaussians  are  rather  poor  in  describing  the 
electronic  behaviour  near  the  nucleus,  since  it  is 
seen  that  a  larger  number  of  functions  are  needed 
to  represent  the  s-orbital  for  comparable  accuracy. 

The  effect  of  a  larger  Gaussian  basis  set  on  the 

computer  requirements  for  molecular  calculations  is 
twofold : 

(a)  The  duration  of  calculation  of  the  electronic  wave  fun¬ 
ction  is  proportional  to  the  number  of  integrals. 

(b)  The  number  of  two-electron  integrals  to  be  evaluated 
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Table  1 .  Convergence  of  the  Hartree-Fock  energy  of  the 


oxygen  atom,  with  the 

number  of  Gaussian 

functions  used  . 

Orbital  No.  of  functions 

Energy  error  (a.u.) 

3 

0.085 

P  4 

0.020 

5 

0.005 

6 

0.001 

7 

0.025 

s  8 

0.008 

9 

0.005 

10 

0.001 

(a)  Reference  11 
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increases  as  the  fourth  power  of  the  number  of 
basis  functions  (i.e.  as  N^)  .  Since  the  time  required 
to  solve  the  matrix  Hartree-Fock  equations  depends  on 
the  number  of  two-electron  integrals,  use  of  large 
Gaussian  basis  sets  can  make  the  iterative  solution 
of  these  equations  very  time  consuming. 

A  number  of  techniques  were  hence  devised  to  ac¬ 
celerate  the  convergence  of  solutions.  These  techniques 
aimed  at  reducing  the  total  time  required,  mainly  by 
reducing  the  number  of  iterations. 

However,  the  most  direct  answer  to  the  time  problem 
should  come  from  an  educated  reduction  of  the  number  of 
basis  functions.  This  was  accomplished  by  contraction, 
that  is  by  grouping  together  some  Gaussian  functions  in  the 
large  basis  set,  and  treating  them  as  one.  Integrals 

are  thus  being  stored  for  this  contracted  function  only. 

12  13 

Successful  contraction  was  shown  "  '  not  to  result  in 
significant  loss  of  energy. 

The  process  of  contraction  is  now  standard  practice. 
The  reduction  of  a  large  Gaussian  basis  set  into  a  smaller, 
compact,  but  flexible  enough  set  of  functions  is  in  the 
heart  of  modern  and  efficient  molecular  calculations. 

It  must  be  stressed  through,  that  a  great  deal  of  caution 
must  be  exercised  when  contracting  atomic  basis  sets. 

There  is  one  more  difference  between  Slater  and 
Gaussian-type  orbitals  that  should  be  mentioned.  The 


•  •  •  9 


functions  behave  as  exp(-^r), 


Slater  Is,  2s,  3s, 

2 

rexp(-Cr),  r  exp(-£r),...  in  contrast,  all  s-type 
Gaussians  are  taken  to  behave  as  exp(-ar^) ,  i.e.  like  1 
functions.  Similarly  all  p-type  functions  behave  like 
the  2p,  etc.  This  was  shown"^  not  to  result  in  any 
adverse  effects  on  the  accuracy  of  calculations. 
Furthermore,  it  simplifies  the  integrals  considerably. 

We  have  so  far  discussed  various  basis  functions, 
but  we  have  not  seen  how  the  atomic  orbital  expansions 
are  made.  This  will  be  the  topic  of  the  next  chapter. 


Ill  ON  STO  EXPANSIONS 


III.l  Gaussian  Expansions  of  Atomic  Orbitals. 

There  have  been  a  number  of  attempts^  to  obtain 
approximate  expansions  of  Slater-type  orbitals  in  terms 
of  Gaussians,  but  it  appears  that  most  of  these  methods 
do  not  produce  satisfactory  results  when  applied  to 
molecular  quantum  mechanics. 

21 

K »  O-ohata,  H.  Taketa  and  S.  Huzinaga  proposed  a 
least  squares  method  of  expansion,  with  special  emphasis 
on  good  results.  A  brief  description  of  their  method, 
which  was  also  used  in  this  work,  is  given  below. 

Define  the  normalized  STO  as: 

’"nXm(^-r)  =  Rn(5;r)  Ylm(e'*)  (  7 

where 


RS(C;r)  =  N^(c)rn  1  exp(-^r) 
n  n 


(  8) 


and 


N®(C)  =  [  (2n)  !  ]  h  (2c)n+!s 
Define  also  the  normalized  GTO  as: 


(  9) 
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(10) 


^nlm (a  ' r>  =  Rn(a;r)  Yim<e'4>> 


where 

R^(a;r)  =  N^(a)rn  1exp(-ar2)  (11) 

and 

Nn(a)  =  2n+1  t  (2n-l)  !  !  ]  ”**  (2t j)~k  a(2n+1>/4  (12) 

The  GTO's  are  restricted  to  have  even  values  of  (n-&-l) . 

In  general  an  N-terra  expansion  of  an  STO  in  terms  of  GTO's 
is  written  as: 


W^r> 


N  G 

-  I  C .  4> .  _  , ,  (a. ;  r ) 

£=1  i  i f  n '  lm  i' 


(13) 


In  equation  (13),  {C^}  are  the  expansion  coefficients. 

The  term  " square  deviation  of  the  Gaussian  expansion", 
is  taken  to  mean  the  deviation  of  the  fitness  of  the 
Gaussian  expansion  with  the  STO;  this  is  measured  by  the 
value  of  the  integral: 


V  (  £ ;  Ci ,  C4^ 


1  -/"■ 


nlm 


N 

(C;r)-E 


i=l 


C  .  (J)? 
iYi  ,n 


lm(ai;r)] 


dV  (14) 


Integration  is  considered  over  the  whole  space. 

Using  the  scaling  theorem,  which  shall  be  treated 
in  a  later  chapter,  the  problem  reduces  to  obtaining  a 
Gaussian  expansion  for  the  STO  with  an  orbital  exponent 
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of  unity. 

The  method  of  least  squares  is  based  on  the  min¬ 
imization  of  the  integral  in  (14),  with  C=l.  This  is 
achieved  by  choosing  a  set  of  {<j>^}  .  The  optimum  set  of 
expansion  coefficients  {C^}  is  determined  by  solving  the 

resulting  linear  equations.  Then,  the  Gaussian  orbital 

.  .  .  .27 

exponents  are  established  using  the  Powell  minimization 
One  of  the  main  difficulties  met  by  this  method  is  the 
existence  of  multiple  minima.  Moreover,  it  is  known  that 
the  number  of  multiple  minima  increases  significantly 
with  the  increasing  number  of  variational  parameters. 
However,  the  results  of  application  of  this  method  are 
judged  good  enough  for  use  in  molecular  calculations. 

In  the  following  section,  one  of  the  best 
established  least  squares  representations  of  Slater-type 
orbitals  as  a  sum  of  Gaussian-type  orbitals  will  be 
discussed . 

III. 2  The  STO-KG  Method. 

The  term  "minimal  basis  set",  is  used  to  indicate  that 
only  one  basis  function  is  used  to  represent  each  occupied 
or  partly  populated  shell  of  the  ground  state  of  an  atom. 
Unpopulated  shells  are  not  represented.  The  contents  of 
the  present  thesis  exclusively  concern  discussions  on 
such  minimal  basis  sets. 

As  already  mentioned,  a  number  of  attempts  have  already 


/■ 


been  made  at  obtaining  GTO  expansions  of  STO's.  Out¬ 
standing  among  these  attempts  is  the  STO-KG  method 

J  O  —  0 

proposed  by  J.A.  Pople  et  al.  .  In  this  method,  each 

atomic  orbital  is  represented  by  an  STO  which  is  least 
squares  expanded  in  terms  of  K  Gaussians.  The  principles 
of  the  method  are  outlined  in  the  following  equations: 


=  £3/2<f/  (l,£r) 


,  K 

4>  (l,r)  =  I  d  .  g,  (a  .,r) 

ns  i_i  ns ,  i  is  ni 


,  K 

4>  (1 , r )  =  l  d  .  g0  (a  .  ,r) 
np  ns,i  2p  ni 


etc . 


(15) 

(16) 

(17) 


There  are  a  number  of  important  features  in  this  scheme, 
(a)  Sis'  g2p'  etc*'  are  the  normalized  Gaussian-type 
functions 


gls (a ’r)  = 

(2cx/tt)  3//4  exp(-ar2) 

(18) 

g2p(a'r)  = 

(128a5/Tr3) 1//4  r  exp(-ar2)  cos  0 

(19) 

etc . 


(b)  The  constants  d  and  a  are  chosen  so  as  to  give  the 
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optimum  least  squares  fit  to  the  corresponding 
STO  with  C=l. 

(c)  The  Gaussian  orbital  exponents  {an^}  are  common  to 
the  ns  and  np  functions.  This  is  done  for  greater 
computational  economy. 

(d)  The  2s,  3s,  etc.,  functions  are  replaced  by  a  linear 
combination  of  Is  functions. 

(e)  All  expansions  correspond  to  an  STO  with  exponent 
C-1.0.  These  exponents  are  later  optimized,  and 
the  scaling  theorem  is  applied  to  the  expansions. 

The  above  method  has  been  extensively  applied  to  a 

number  of  atomic  and  molecular  systems.  It  enjoys  a 
fair  amount  of  success  in  predicting  the  geometrical 
structures  of  many  neutral  molecules  in  their  ground 
states . 

Increasing  the  value  of  K,  i.e.  increasing  the  size 
of  the  expansion,  makes  the  results  of  the  STO-KG  method 
approach  those  of  the  Slater-type  basis  itself.  It  is 
obvious,  however,  that  the  method  can  not  surpass  the 
accuracy  of  STO-SZ  (STO  sincrle-zeta)  calculations.  Thus, 
it  is  unfortunate  that  STO-KG  is  expected  to  be  inappropriate 
for  use  in  treating  third  row  atoms  due  to  this  inherent 
limitation.  This  point  will  be  raised  again  and  discussed 
to  a  greater  depth  in  a  later  chapter. 

The  next  chapter  introduces  a  new  type  of  minimal 
Gaussian  basis  set,  which  aims  at  overcoming  some  of  the 
limitations  of  the  STO-KG  method. 


/■ 
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IV  A  NEW  MINIMAL  GAUSSIAN  BASIS  SET 


The  minimal  Gaussian  basis  set  proposed  in  this  work 
aims  at  better  orbital  energies,  with  particular  emphasis 
on  the  bonding  valence-orbitals  of  atoms  for  more  accurate 
quantum  mechanical  molecular  calculations.  The  orbital 
expansions  still  remain  at  the  three-term  (3G)  level,  but 
a  new  approach  to  their  formation  is  suggested. 

IV. 1  The  Atomic  Orbital  Expansions. 

Each  atomic  orbital  will  be  expressed  as  a  linear 
combination  of  three  Gaussian-type  functions.  The 
expansion  method  is  the  following: 

For  the  Is  orbital,  a  three-term  expansion  is  written  as: 


3 


(20) 


2  8 

This  expansion  is  fitted  against  the  reference  STO 

double-zeta  Is  function,  according  to  the  least  squares 

21 

fitting  scheme  proposed  by  S.  Huzinaga  et  al .  and 
discussed  in  section  III.l.  The  optimum  values  of  the 
constants  an  .  and  C,  ,  are  thus  determined. 


Is ,  K 


In  a  completely  analogous  manner,  we  expand  the  2s 


orbital  as: 


3 


(21) 


20 


' 
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To  the  expansion  (21)  we  now  incorporate  the  contracted 

function  from  (20)  .  This  augmentation  of  (21)  has  a  very 
important  purpose.  Inclusion  of  the  <J>^g  serves  to 
describe  the  inner  part  of  the  2s  orbital,  so  that  the 
a0  ,  and  C0c,  ,  are  freer  to  describe  better  the  outer 
part  of  the  orbital.  The  2s  expansion  now  becomes: 


2s 


=  C 


<|)ls 


3 

^  ^2  s  k 

k=l  ' 


g2  s ,  k  ^a2s ,k' 


(22) 


It  is  this  expansion  (22)  that  we  now  fit  against  the 
corresponding  STO  double-zeta  function  for  the  atom. 

Thus,  the  optimum  values  of  the  constants  o^g  k  and  C 2S ^ 
are  determined. 

2  8 

We  now  refit  the  Is  expansion  to  the  STO  double-zeta 
Is  function  for  the  atom.  This  time,  however,  we  include 

in  the  expansion  the  contracted  cf> ^ s  from  (21)  .  As  before, 

the  contracted  2s  will  describe  the  outer  part  of  the  Is 

orbital,  so  that  the  als  k  and  Clg  k  are  freer  to  describe 

better  the  inner  part.  The  new  Is  orbital  expansion  now  is: 


4> 


Is 


3 

Z  C,  . 
k=l  1S'k 


gls,k (als,k'r) 


(23) 


From  the  fitting  we  obtain  a  new  set  of  (ot^  k*  *^ls,k^  f 
k=l-3 . 

The  expansions  (22)  and  (23)  are  now  fitted  suc¬ 
cessively  against  the  STO  double-zeta  corresponding 
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functions  ,  until  the  resulting  least  squares  standard 
deviations  are  minimized. 

The  2p  orbital  is  similarly  expanded  and  fitted. 


4> 


2p 


3 

E  ,  g~  ,  (a0  ,  ;r) 

k=1  2p,k  y2p,k  2p,k' 


(24) 


If  3p  or  higher  orbitals  are  present,  their  method  of 
treatment  is  completely  analogous  to  that  described  for 
the  Is  and  2s  functions. 

At  the  end  of  the  fitting  process,  there  exists  a 

set  of  three  (a  ,  ,  C  ,  )  pairs,  that  completely  charac- 

n ,  x  n ,  x 

terize  the  three-term  expansion  for  each  atomic  orbital. 
1.00 


Is : 

[ (als,l;Cls,l)  ' 

(als,2;Cls,2) ' 

(als,3;Cls,3) ] 

2s  j 

[ {a2s,l;C2s,l)  ' 

(a2s , 2 ' ^2s , 2  ^ ' 

^a2s , 3 ' *^2s , 3  ^ ^ 

2p : 

[ (a2p,l;C2p,l) ' 

(a2p,2;C2p,2) ' 

(a2p,3;C2p,3) J 

etc  o 

Each  of  the  above  sets  provides  a  contracted  basis 
function  which  will  be  further  energetically  optimized 
by  the  method  described  in  the  next  section. 


IVc2  Energetic  Optimization  of  the  Orbital  Expansions. 

The  atomic  programme  used  for  all  atomic  calculations 
is  a  standard  Hartree-Fock-Roothaan  SCF  programme,  written 
in  Fortran  by  S.  Huzinaga,  and  based  on  the  Roothaan 
expansion  method. 

To  optimize  the  basis  set  as  determined  in  the 
previous  section  and  reported  in  the  sets  (25) ,  a  standard 
method  was  employed  for  consistency.  We  shall  describe 
this  method  for  atoms  having  Is,  2s,  and  2p  (fully  or 
partially  occupied)  shells.  The  procedure  consists  of 
the  following  steps: 

Stel  1:  Perform  an  atomic  calculation,  with  all  basis 
functions  contracted.  The  basis  set  in  this 
case  is:  GTO(33/3). 

Step  2:  Release  the  contracted  function  for  the  Is 
orbital,  keeping  all  the  others  contracted. 

Thus,  optimize  the  Is  orbital  expansion  co¬ 
efficients  and  orbital  exponents.  The  basis 
set  in  this  case  is:  GTO(1113/3). 

Step  3:  Using  the  orbital  exponents  and  contraction 

coefficients  obtained  in  step  2,  contract  the 
Is  orbital  and  release  the  2s  to  optimize  the 
a0  .  and  C0  ,  ,  k=l-3.  In  this  case,  the 
basis  set  is:  GTO(3111/3). 

Recontract  the  2s  orbital  with  the  contraction 


Step  4 : 
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coefficients  and  orbital  exponents  determined 
in  step  3.  Release  to  optimize  the  parameters 
of  the  2p  orbital.  The  basis  set  is:  (33/111). 
Step  5:  Using  the  values  of  a ^  ^  and  k  0^ta^-ne(^  an 

step  4,  recontract  the  2p  orbital.  With  all 
basis  functions  contracted,  perform  an  atomic 
calculation.  This  step  is  the  same  as  step  1, 
and  the  basis  set  is  GTO(33/3). 

If  3s  and  3p  orbitals  also  exist,  insert  analogous 
steps  after  steps  3  and  4  respectively.  If  in  addition, 
the  3rd  orbital  is  also  partly  or  wholly  occupied,  a 
similar  step  is  inserted  before  step  5.  The  same 
reasoning  applies  for  even  higher  orbitals. 

At  the  next  stage  of  the  process,  steps  2  through  to 
5  are  successively  repeated  in  cycles.  At  the  end  of 
each  cycle  the  value  of  the  atomic  energy  drops,  as 
compared  to  the  previous  cycle.  The  repetition  of 
these  cycles  stops,  when  the  drop  in  energy  from  one 
cycle  to  the  next  becomes  insignificantly  small.  It 
was  found  that  three  to  four  such  cycles  were'  sufficient 
for  calculations  with  the  first  row  atoms.  For  atoms  of 
the  second  and  third  periods,  two  to  three  cycles  were 
judged  enough  to  provide  the  desired  results. 


The  values  for  the  orbital  exponents  and  contraction 
coefficients  determined  by  the  above  method  are  given  in 
Table  2  for  the  first  row  atoms  in  their  ground  states. 


' 


Number  in  parenthesis  denotes  power  of  10. 
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Table  2;  Orbital  Exponents  and  Coefficients  for  the  first  row  atoms. 


In  Table  3  the  values  of  these  parameters  are  given  for 
Sulphur  and  Nickel. 

IV. 3  Atomic  Calculations:  Results  and  Discussion. 

In  the  preceding  section  we  arrived  at  a  complete 
basis  set  for  each  atom  of  the  first  row  of  the  periodic 
table,  and  also  for  Sulphur  and  Nickel  at  their  ground 
states.  These  basis  sets  have  been  given  in  Tables  2 
and  3.  The  atomic  orbital  energies  determined  using  the 
above  sets  are  displayed  at  Tables  4  and  5.  In  these 
tables,  the  values  of  STO-SZ,  STO-DZ  and  H.F.  calculations 
are  also  given  for  reasons  of  comparison. 

We  have  also  calculated  the  percent  relative  error 
of  our  atomic  orbital  energies  with  the  corresponding 
STO-SZ  and  H.F.  results.  These  are  given  by 


respectively.  Their  values  are  displayed  in  Table  6.  A 
minus  (-)  sign  before  the  number  means  that  our  calculated 
orbital  energy  is  lower  than  the  corresponding  STO-SZ, 
by  the  percentage  indicated  by  that  number . 

A  number  of  comments  may  be  made  on  observation  of 
Table  6 . 

(a)  All  calculated  Is  orbital  energies  are  poorer  than 


Table  3:  Orbital  Exponents  and  Coefficients  for  S  and  Ni . 


Orb 

Sulphur 

(3P) 

Nickel 

(3f) 

Exponent 

Coefficient 

Exponent 

Coefficient 

1.12075 (3) 

7.00404 (-2) 

3.71633(3) 

6.46093 (-2) 

Is 

1.72613(2) 

3.90016 (-1) 

5.67873(2) 

3.74593 (-1 ) 

3.84041(1) 

6 .51111 (-1) 

1.25777(2) 

6.64236 (-1) 

4.84113(1) 

-9.99127 (-2) 

1 .56865 (2) 

-1 .11233 (-1) 

2s 

5.23960 

4. 87790 (-1) 

1.82921(1) 

5.37274 (-1) 

2.15444 

6.06343(-l) 

8.12609 

5 .55694 (-1) 

2.52106 

"1.76859 (-1) 

1.36957(1) 

-2.32084 (-1) 

3s 

4.03651 (-1) 

6. 77064 (-1) 

2.25807 

5 .77746 (-1) 

1.47347 (-1) 

5  .01158  (-1) 

1.02896 

6.29956 (-1) 

1.36146 

-1 . 87838 (-1) 

4s 

1.24272 (-1) 

5.62553(-l) 

4 .83679 (-2) 

5.45291 (-1) 

4.54572(1) 

1 .28901 (-1) 

1.98521(2) 

1 . 00796 (-1 ) 

2p 

1.04817(1) 

4 .86214 (-1) 

4.64787(1) 

4 . 56571 (-1) 

2.91319 

5 .53493 (-1) 

1.36427(1) 

0.77214 (-1) 

7.53257 (-1) 

3.04283 (-1) 

3.43541(1) 

-2.43491 (-2) 

3p 

2.89607 (-1) 

5. 10095  (-1) 

3.90276 

5 . 15717  (-1) 

1.10548 (-1) 

3.  46945 (-1) 

1.20099 

6 . 38098 (-1) 

1.32456(1) 

1.89371 (-1) 

3d 

3.34272 

5.22981 (-1) 

8.10345 (-1) 

5  .66377 (-1) 

Number  in  parenthesis  denotes  power  of  10. 
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Table  4:  Orbital  Energies  (  a.u.)  for  the  first  rov:  c*toms 
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Table  5:  Orbital  Energies  (a.u.)  for  Sulphur  and  Nickel 


Orb 

Calc. 

s(3r) 

Ni(3F) 

Present 

-91.40960 

-303.62705 

Is 

STO-S Za 

-91.08467 

-305.50899 

STO-DZ3 

-92.00112 

-305.58756 

„  „  a 

H .  F . 

-92.00447 

-305.61S82 

Present 

-8.98057 

-37.72615 

2s 

STO-SZ3 

-8.77098 

-37.21212 

STO-DZ3 

-9.00305 

-37.89093 

„  „  a 

H.F. 

-9.00431 

-37.91793 

Present 

-0.87G56 

-4.79509 

3s 

STO-S Z3 

-0.81168 

-4.34110 

STO-DZa 

-0.87897 

-4.84609 

_  a 

H  •  F  • 

-0.87955 

-4.88783 

Present 

-0.25482 

4s 

STO-SZ3 

-0.22276 

STO-DZ3 

-0.26640 

H.F. 3 

-0.27624 

Present 

-6.62400 

-32.71390 

2p 

STO-SZ3 

-6.43647 

-32.63106 

STO-DZ3 

-6.68121 

-32.89703 

a 

HcF. 

-6.68253 

-32.94182 

Present 

-0.43356 

-3.17493 

3p 

STO-SZ3 

-0.38720 

-2.85015 

STO-DZ3 

-0.43694 

-3.23877 

H.F.3 

-0.43737 

-3.27764 

Pre  sent 

-0.54799 

3d 

STO-SZ3 

-0.00054 

STO-DZ3 

-0.66185 

H.F.3 

-0.70693 

(a)  Ref.  28 
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STO-SZ  by  an  average  of  0.58%. 

(b)  All  the  other  orbital  energies,  (excluding  the  2s  of 
Be),  are  superior  to  STO-SZ. 

(c)  For  the  first  row  atoms,  the  accuracy  of  the  2s 
energy  value  increases  with  atomic  number  as  compared 
to  STO-SZ.  Beryllium  is  again  excluded.  The 
percent  superiority  over  STO-SZ  varies  from  -0.29% 
for  Lithium,  to  -9.11%  for  Neon. 

(d)  The  same  trend  is  obvious  for  the  2p  orbital  energies, 
where  the  percentage  improvement  varies  from  -0.69% 
for  Boron,  to  -40.3%  for  Neon. 

(3)  For  the  first  row  atoms  again,  it  is  noted  that  the 
accuracy  of  the  orbital  energies  as  compared  to 
STO-SZ  increases  as  one  proceeds  to  higher  orbitals. 
Consider  the  example  of  Fluorine.  The  Is  calculated 
value  is  0.54%  poorer  than  STO-SZ;  the  2s  is  7.92% 
better,  and  the  2p  is  30.1%  better  than  STO-SZ. 

This  is  an  extremely  important  feature  because  one 
is  more  concerned  about  the  accuracy  of  the  outer 
bonding  orbitals,  since  these  are  the  ones  that  mainly 
participate  in  molecular  bond  formation. 

(f)  Similar  trends  are  observed  for  Sulphur  and  Nickel. 

The  most  important  feature  presented  here  is  the 
value  of  the  Nickel  3d  orbital  energy.  Whereas  STO-SZ 
gives  a  value  of  -0.00054  a.u.  and  ST0-3G  a  value  of 
+0.05211  a.u.,  our  calculated  value  is  -0.54799  a.u. 


to  be  compared  with  the  Hartree-Fock  result  of 
-0.70693  a.u.  Further  comments  on  this  important 
point  will  be  made  later  on  in  this  section. 

The  atomic  energies  obtained  using  our  basis  set  are 
displayed  in  Table  7.  In  the  same  table  the  atomic 
energies  of  ST0-3G,  STO-4G,  STO-SZ,  STO-DZ  and  H.F.  are 
given  for  comparison.  With  respect  to  these  values,  the 
following  comments  can  be  made. 

(a)  Our  calculated  atomic  energies  are  invariably  higher 
than  STO-SZ.  However,  the  difference  appears  to  be 
fairly  small. 

(b)  All  our  values  are  lower  than  STO-3G.  The  dif¬ 
ference  keeps  growing  to  the  favor  of  the  present 
method,  with  growing  atomic  number. 

(c)  For  Oxygen  the  atomic  energy  calculated  with  the 
basis  proposed  in  this  work  approaches  the  value  of 
STO-4G .  For  all  heavier  atoms  our  method  yields 
results  that  are  superior  to  those  of  STO-4G. 

(d)  For  Nickel,  the  atomic  energy  calculated  here  is 
deeper  than  that  of  STO-4G  by  approximately  5.72 
hartrees . 

In  Table  8  we  calculate  the  relative  percent  error 
of  our  atomic  energy  results  as  compared  to  those  of 
STO-SZ  and  H.F.  The  formulae  used  are: 


respectively . 


' 


Table  7:  Atomic  Energies  (a.u.)  of  the  first  row  atoms,  Sulphur  and  Nickel 
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Table  8.  Relative  percent  error  of  calculated  atomic 

energies  as  compared  to  STO-SZ  and  H.F.  values. 


Atom 

STO-SZ 

H.F. 

Li 

0.55 

0.74 

Be 

0.75 

0.69 

B 

0.53 

0.65 

C 

0.46 

0.64 

N 

0.39 

0.63 

0 

0.28 

0.63 

F 

0.17 

0.64 

Ne 

0.08 

0.65 

S 

0.29 

0.51 

Ni 

0.12 

0.44 

(a)  Ref.  28 


The  average  percent  error  on  comparison  with  the 
Hartree-Fock  values  is  0.62%.  The  pattern  of  the  errors 
resulting  from  the  comparison  with  the  STO-SZ  values 
indicates  that  the  accuracy  of  our  method  improves  with 
increasing  atomic  number.  Again  the  results  for  Nickel 
deserve  particular  attention. 

Table  9  displays  a  comparison  between  STO-SZ  and 
Hartree-Fock  energies  for  the  3d  orbitals  of  the  third 
row  atoms.  It  is  obvious  that  STO-SZ  is  extremely  in¬ 
effective  in  treating  these  atoms.  The  case  of  Copper 
clearly  demonstrates  the  bad  and  unrealistic  results 
obtained  by  the  use  of  STO-SZ  in  these  cases.  Naturally, 
one  would  not  expect  STO-3G  to  perform  better  than  the 
single-zeta  method.  The  unavoidable  conclusion,  therefore, 
is  that  STO-3G  should  not  be  used  in  calculations  in¬ 
volving  transition  metals. 

In  summary,  we  can  say  that  our  method  invariably 
gives  better  atomic  energies  than  STO-3G.  In  addition, 
our  orbital  energies  are  also  better,  with  particular 
emphasis  on  the  accuracy  of  the  higher  out-lying  atomic 
orbitals.  We  would  thus  expect  to  obtain  better  molecular 
results.  Most  important,  the  poor  performance  of  STO-3G 
and  STO-SZ  in  treating  transition  metals  has  been  proven 
beyond  any  question.  Hence,  it  can  safely  be  concluded 
that  our  method  is  highly  appropriate  for  calculations 
involving  third  row  atoms. 


. 
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Table  9:  STO-SZ  and  H.F.  3d  orbital  energies  (a.u.)  for 
3rd  row  atoms. 


STO-SZa 

H.F.a 

.2  % 

Sc  (  D) 

-0.18988 

-0.34357 

Ti (3F) 

-0.24338 

-0.44071 

v(4f) 

-0.25642 

-0.50966 

Cr  (  D) 

-0.25083 

-0.56920 

Mn  (6S^ 

-0.25159 

-0.63884 

,5  . 

Fe  (  D) 

-0.15690 

-0.64689 

,4  x 

Co(  F) 

-0.08085 

-0.67548 

:Jl(3f) 

-0.00054 

-0.70693 

,2  x 

Cu(  S) 

0.50637 

-0.49074 

(a)  Ref.  28 


V  MOLECULAR  CALCULATIONS 


V.l  The  Molecular  Programme. 

29 

For  all  our  molecular  calculations,  the  MOLECULE 
package  programme  has  been  used.  This  programme  performs 
calculations  within  the  Hartree-Fock  S.C.F.  scheme.  The 
integral  calculations  assume  an  atom  centered  GTO  basis 
set . 

For  our  purposes  the  programme  consists  of  two 
separate  parts: 

(a)  Integral  Calculation  (MOLECULE) .  The  integral 

30 

package  was  written  by  J.  Almlof  ,  but  it  includes 
a  number  of  modifications  due  to  O.  Gropen.  It 
also  includes  the  option  of  using  model-potentials 

31 

to  simulate  the  inner-core  according  to  S.  Huzinaga 
The  programme  efficiently  uses  the  molecular  symmetry 
in  the  evaluation  of  two-electron  integrals.  Thus, 
recalculation  of  symmetry-redundant  integrals  is 
avoided,  and  the  construction  of  blocks  of  symmetry 
matrices  and  supermatrices  is  achieved  without  the 
time-consuming  transformation  procedure. 

(b)  Closed  shell  S.C.F.  This  is  a  standard  S.C.F. 
closed  shell  programme.  It  has  options  for  trial 
vectors,  extrapolation  procedure,  and  use  of  the 
SHIFT-operator  and  damped  F-matrix  procedure  for 
obtaining  a  better  convergence. 
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In  Figures  1  and  2  we  give  the  approximate  con¬ 
structions  of  the  two  parts  of  the  programme-package 
described  above. 

V.2  The  Scaling  Theorem. 

Reference  to  this  theorem  has  already  been  made  in 
chapter  III. 2.  It  will  also  be  used  again  in  the  next 
section.  This  is,  therefore,  an  appropriate  point  to 
give  an  account  of  the  theorem. 

The  scaling  theorem  can  be  stated  as  follows: 

Suppose  that 

„  N 

\p  (l;r)  =  E  C.  <j>.  (a.  ;r)  (26) 

i=l  1  1  1 

is  the  Gaussian  expansion  of  a  Slater-type  orbital,  with 

the  exponent  value  of  unity.  Then,  the  Gaussian  expansion 

for  the  STO  with  an  arbitrary  exponent  value  of  £ ,  will 

have  the  same  expansion  coefficients  {C^} ,  but  all  the 

Gaussian  exponent  parameters  {ou}  will  be  multiplied 

2 

(scaled)  by  a  factor  of  £  .  Thus, 

N  r  9 

\pS  (  £  ;  r )  =  I  ci(j)V(£zai;r)  (27) 

i=l 

We  will  first  give  a  dimensional  analysis  of  the 
theorem.  Let  us  rewrite  equation  (26)  as: 


, 


Figure  1 


Construction  of  the  Integral  Part 


READIN:  Read  input,  create  necessary  symmetry  information 

and  normalize  contracted  basis  functions  and 
symmetry  functions. 

CORPOT:  Read  input  for  core-part  when  model-potential 

option  is  used.  Organize  the  core-part  in¬ 
formation.  Use  a  special  subroutine  for  the 
calculation  of  overlap-integrals. 

ONEL:  Calculate  overlap  kinetic  energy  and  nuclear 

attraction  energy  integrals. 

Calculate  two-electron  integrals. 


TWOEL: 


' 


Figure  2 


Construction  of  the  S.C.F.  Part 


SCFIN :  Read  input  data,  organize  and  print. 

ITERAT:  Organize  the  S.C.F.  procedure. 
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(28) 


_  N 

^  (r )  =  Z  C .  (K  (a  .  r  ) 
i=l  1  1  1 


Here,  we  have  made  the  STO  and  GTO  functions  depend  on 

one  parameter  only  instead  of  two.  This  is  perfectly 

legitimate  since  the  exponential  part  of  the  STO  with 

C=1  is  exp(-r);  similarly  the  exponential  part  of  the 

2 

GTO  is  exp(-ar  ). 

2 

We  now  make  the  mapping  cu  **■  X>  cu  ,  where  c  is  a 
real  number.  Then: 


Z  C .  cf>9  ( a.r2)  Z  C .  <j>?  ( £2a  .  r2 )  = 
i=l  11  1  i=l  11 


N  ? 

Z  C.<|>.  [a.  Ur)  ] 
.  ,  1Y1  1 
1=1 


(29) 


By  observing  equation  (28) ,  we  can  say  that  (29) 
corresponds  to  the  left-hand  side  of  (28) ,  if  the  latter 
is  written  as: 


( tr) 


-1  2 

By  necessity,  C  should  have  dimensions  of  L  ,  and  £ 


should  be  L 


-2 


A  proof  of  the  theorem  is  given  in  Reference  21.  We 


shall  consider  here  a  particular  feature  that  makes  the 
transition  from  (26)  to  (27)  possible.  Throughout  this 


42 


discussion  it  must  be  kept  in  mind  that  the  coefficients 

2 1 

(C^)  have  been  shown  to  be  independent  of  C 

Consider  the  right-hand  side  of  relation  (26)  and 

take  the  overlap  integral.  Since  we  are  dealing  with 

orthonormalized  Gaussian  functions,  one  only  has  to 

examine  the  cross-terms  <§./§.>.  The  definitions  of 

13 

G  C 

<|k  and  <Jk  are  taken  from  equation  (10)  in  chapter  III.1. 
Then : 


00 


13 


=  [N^(a) ]2  J  r2 (n  1}e  (ai 


+a . ) r  2  , 

1  j  r  dr 


22 (n+i)  a(2n+l)/2  [ (2n~l) ! I ] 

- r 


TT 


T7l  (a- a.) 


(2n+l)/4 


[  ( 2n-l )  1  !  ]  (2TT)'*  2 

^  (n+h) 


0  (n+1)  '10 


(a . +a . )  V  (a . +a . ) 

1  3  V  1  0 


2 (a^a j ) 


a  .  +  a  . 
1  3 


(30) 


Equation  (30)  is  worth  commenting  upon.  Suppose  that 
instead  of  cu  and  a  ,  we  had  £  and  6  a  j  •  Then  it  is 
obvious  that  the  factor  6  2  would  be  eliminated  from  (30) . 
This  fact,  combined  with  the  fact  that  the  coefficients 
(Ci)  are  independent  of  C,  make  the  scaling  theorem 
valid . 


V.3  The  Method. 

The  basis  sets  of  Tables  2  and  3  were  used  in  a 


number  of  calculations  on  molecular  systems.  It  was 
found  that  although  these  calculations  yielded  reasonable 
values  of  molecular  and  molecular  orbital  energies,  they 
did  not  produce  any  binding  energies.  This  defect  was 
attributed  to  the  fact  that  the  basis  sets  developed  by 
atomic  calculations,  failed  to  take  into  account  the 
effects  of  molecular  formation  on  the  atomic  orbitals. 

The  need,  therefore,  arose  of  scaling  the  Gaussian 
functions  used  in  these  basis  sets.  For  this  purpose, 
the  scaling  theorem  discussed  in  the  preceding  section  was 
used . 

With  the  aim  of  keeping  the  method  at  a  fairly 
simple  level,  it  was  decided  to  scale  only  the  outer 
valence  atomic  orbitals.  Calculational  experimentation 
indicated  that  this  practice  did  not  seem  to  significantly 
worsen  the  results.  In  accordance,  only  the  2p  orbital 
was  scaled  for  the  first  row  atoms,  and  the  3p  orbital 
for  Sulphur.  The  optimum  scaling  factor  of  an  atom  in  a 
particular  molecule  was  determined  by  successive  calculat¬ 
ions  until  the  total  molecular  energy  was  a  minimum. 

From  a  series  of  such  optimizations  of  different  molecules, 
an  average  scale  factor  for  the  valence  orbital  of  that 
atom  was  established.  The  average  optimum  scaling 
factors  for  hydrogen,  carbon,  nitrogen,  oxygen  and 
sulphur  are  reported  in  Table  10. 

All  molecular  calculations  were  performed  at  the 


.  J 
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Table  10:  Valence  orbital  average  scaling  factors 


Atom 

Orb 

- -  ....  - - --■■■■  ...  .  .  ■— ■  — > 

Scale  factor 

H 

Is 

1.450 

C 

2p 

1.065 

N 

2p 

1.065 

0 

2p 

1.065 

S 

3p 

1.100 

-  ■■  -  - - - J - — 

, 

geometries  shown  in  Table  11.  In  most  cases  these 
structural  parameters  are  Hartree-Fock  equilibrium 
geometries . 

33 

In  parallel, the  package  programme  GAUSSIAN-70 

was  used  to  calculate  the  same  molecules,  with  an  ST0-3G 

basis  set,  at  the  above  values  of  the  geometrical 

parameters  (Table  11) .  This  was  done  in  order  to  be 

able  to  compare  the  results  of  the  two  methods  at  the 

same  geometry  for  a  certain  molecule.  It  must  be  noted 

that  when  a  geometry  optimization  was  performed  with 

GAUSSIAN-70  and  ST0-3G,  the  results  were  the  same  as  those 

„  22-26 

reported  by  J.A.  Pople  et  al . 

We  also  carried  out  a  geometry  optimization  of  the 
molecular  systems  we  treated.  For  diatomic  molecules 
this  was  done  by  the  least  squares  fitting  of  the  energy 
curve  via  a  third— degree  polynomial.  For  triatomic  or 
polyatomic  molecules,  the  geometrical  parameters  were 
optimized  step  by  step  and  one  by  one,  making  use  of  the 
symmetry  point  group,  to  which  the  molecule  belongs. 

In  the  next  section,  we  shall  present  the  results  of 
our  molecular  calculations.  A  discussion  of  these 
results  and  a  comparison  with  other  methods  will  also  be 


given . 


* 


/ 
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Table  11:  Molecular  Geometries  for  Table  12 


Molecule 

r (a. u  . ) 

0 (deg . ) 

N2 

r  =  2.0G80 

C2 

r  =  2.3481 

CO 

r  =  2.1320 

0-2 

r  =  2.1944 

0  =  180.0 

H2° 

r  =  1.8085 

0  =  104.5 

S°2 

r  =  2.7061 

0  =  119.5° 

HCN 

r  =  2.0143 

C-H 

r  =  2.1791 

C-N 

e  =  180.0 

C2»2- 

r  =  2.2810 

c-c 

r  =  2.0020 

C-H 

0  =  180.0 

HCC 

47 


V.4  Results  and  Discussions. 

The  average  values  of  the  atomic  valence  orbital 
scaling  factors  appear  in  Table  10.  Table  11  contains 
the  geometries  of  the  ground  states  of  the  molecules  used 
in  the  calculations.  All  results  given  in  Table  12 
correspond  to  these  geometrical  parameters. 

The  STO-3G  results  also  displayed  in  Table  12,  were 

33 

calculated  here,  using  GAUSSIAN-70  .  The  values 
reported  are  at  the  geometries  of  Table  10. 

It  can  be  seen  that  our  calculated  molecular 
energies  are  consistently  lower  than  those  of  STO-3G. 

For  sulphur  dioxide  the  difference  between  the  two  methods 
is  approximately  3.2  hartrees.  In  molecules  involving 
first  row  atoms,  our  method  gives  energies  lower  by  an 
average  of  0.62  hartrees,  than  STO-3G.  On  comparison 
with  the  Hartree-Fock  or  near  Hartree-Fock  molecular 
energies,  the  average  relative  error  of  our  method  comes 
to  0.75%  for  the  molecules  tested.  Upon  the  same 
comparison,  the  STO— 3G  produces  results  with  an  average 
relative  error  of  1.34%. 

In  Table  12,  the  binding  energies  are  also  given 
underneath  the  molecular  energies.  In  most  cases  they  are 
of  the  same  order  of  magnitude  as  those  of  STO-3G. 
Exceptions  are  the  N2  end  CC>2  molecules,  where  they  appear 
to  be  fairly  small  in  absolute  value.  For  C2  and  SC^ 


i  g  *  .j-  .  <  ■ 
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Table  12: 

Molecular  and 

Binding  Energies 

(a .u . ) 

Molecule 

This  Method 

STO-3G6 

nr.  H.F. 

No 

-108.12076 

-107.49484 

-108.99283 

2 

-0.00506 

-0.05464 

-0.1909 

c_ 

-74.85262 

-74 .42205 

-75.40b20a 

2 

+0.04424 

+0.03527 

-0.02898 

CO 

-111.90495 

-111.22458 

-112.77901b 

-0.12248 

-0.19167 

-0.28103 

C0o 

-186.19402 

-185.06472 

-187.70366b 

2 

-0.07751 

-0.22756 

-0. 39631 

HO 

-75.46393 

-74.96291 

-76 . 05936° 

2 

-0.13593 

-0.16884 

-0.25000 

S0o 

-543.79626 

-540.60254 

-547.22464d 

2 

+0.36199 

+0.18491 

-0.12545 

HCN 

-92.19859 

-91.67520 

b 

-92.90868 

-0.19533 

-0.23153 

-0.31914 

C  II 

-76.28170 

-75.85203 

-76.84959b 

2  2 

-0.39088 

-0.40489 

-0.47237 

(a) 

Ref. 

34 

(b) 

Ref. 

35 

(c) 

Ref . 

37 

(d) 

Ref. 

36.  Also  HF(DZ)  values:  E 

-546.99113,  Ebind  =  +0.09364. 


(e)  Our  own  calculation  using  GAUSSIAN- 70,  with  STO-3G  basis  set. 
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neither  our  method  nor  ST0-3G  produce  any  binding  energy. 

3  6 

It  should  be  stressed,  however,  that  the  H.F.  (DZ)  method 
also  yields  no  binding  for  sulphur  dioxide. 

The  results  of  the  geometry  optimization  are  given 
in  Table  13,  which  displays  the  calculated  equilibrium 
geometries.  With  the  exception  of  C02 ,  all  other 
structural  parameters  are  closer  to  the  near  Hartree- 
Fock  or  experimental  values  than  ST0-3G. 

Finally,  in  Table  14  we  give  the  molecular  orbital 
energies  of  N2  and  C02  at  the  internuclear  distances  of 
Table  11.  The  ST0-3G  values  have  been  obtained  using 
GAUSSIAN-70  at  the  above  geometries.  Our  results  are 
again  closer  to  the  Hartree-Fock  energies  than  ST0-3G. 

The  difference  between  the  molecular  orbital  energy  values 
yielded  by  the  two  methods  is  distinctly  pronounced  for 
the  outer  orbitals.  This  is  also  to  the  advantage  of  the 
present  method,  since  most  of  the  molecular  properties 
depend  on  the  location  of  these  outer  orbitals. 


Table  13:  Calculated  Molecular  Geometries;  r(a.u.),  G(deg.). 


.  Molecule 

This  method 

STO- 3G 

(nr) H .F. 

Expt . 

N2 

r  =  2.1383 

2.3300d 

2.0680a 

2.0674d 

C 

r  =  2.4518 

2.1429d 

2.34813 

2.3470d 

2 

CO 

r  =  2.1112 

2.1645a 

2.1320b 

2.1320 

C°2 

r  =  2.2900 

2.2450d 

2.1944b 

2.1921d 

6  =  180.0 

180.0 

180.0 

180.0 

H  O 

r  =  1.8347 

1.8708d 

1.800e 

1.8085d 

2 

0  =  102.5 

100.0 

105.0 

104.5 

s°2 

r  =  2.7300 

2.8911f 

2.9329a 

2.7061° 

0  =  112.5 

119.45 

112.4 

119.5 

d 

b 

_  d 

HCN 

r  =  2.0163 

2.0220 

2.0143 

2.0088 

C-H 

r  =  2.1792 

2.1788 

2.1791 

2.1807 

C-N 

0  =  180.0 

180.0 

180.0 

180.0 

C  H 

r  =  2.2604 

2.2072d 

2.2810b 

2.1733d 

2  2 

C-C 

r  =  2.0005 

2.0126 

2.0020 

2.0050 

C-H 

0  =  180.0 

180.0 

180.0 

180.0 

HCC 

(a)  Ref.  34 

(b)  Ref.  35 

(c)  Ref.  36 

(d)  Ref.  24 

(e)  Ref.  37 

(f)  Our  own  calculation  using  GAUSSIAN-70  and  STO-3G  basis  set 
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Table  14:  Molecular  Orbital  energies  (a.u.)  for  N2  and  CO 


N2  (r 

=  2.0680  a 

.u. ) 

CO  (r 

=  2.1320  a. 

u.) 

Orb 

a 

STO-3Gb 

H.F.C 

Orb 

a 

b 

STO-3G 

H.F.d 

la 

g 

-15.58252 

-15.51928 

-15.63195 

la 

-20.44436 

-20.42433 

-20.67389 

la 

u 

-15.58389 

-15.51727 

-15.67833 

2a 

-11. 32088 

-11.09336 

-11.36987 

2  a 

g 

-1.48374 

-1.44619 

-1.47360 

3a 

-1.52434 

-1.46006 

-1.52826 

2a 

u 

-0.75550 

-0.72202 

-0.77796 

4a 

-0.76506 

-0.69951 

-0. 81091 

ITT 

U 

-0.62031 

-0.57554 

-0.61544 

5a 

-0.51902 

-0.44645 

-0.55970 

3a 

g 

-0.58456 

-0.54039 

-0.63495 

6  a 

0.94367 

1.04021 

0.42782 

ITT 

g 

0.23248 

0.282S0 

ITT 

-0.62939 

-0.55117 

-0.64529 

3a 

u 

1.15783 

1.13178 

2ir 

0.21260 

0.31344 

0.16542 

(a)  Present  Method 

(b)  Our  own  calculation  with  GAUSSIAN-70  and  ST0-3G  basis  set. 

(c)  Ref.  34 


(d)  Ref.  35 


VI  CONCLUSION 


The  objective  of  the  present  work  was  to  develop  a 
new  minimal  Gaussian  basis  set  to  be  used  in  atomic  and 
molecular  calculations.  The  basis  set  had  to  be  at  a 
fairly  simple  level,  so  that  it  could  be  versatile  and 
useful.  Since  computer  time  is  a  major  factor  to  be 
taken  into  consideration  with  most  modern  quantum 
mechanical  calculations,  the  basis  set  had  to  be  ec¬ 
onomical  in  terms  of  time  when  applied  to  such  calculations. 
The  atomic  orbital  expansions,  therefore,  had  to  be  made 
in  terms  of  a  reasonably  small  number  of  Gaussian 
functions.  This  is  a  common  decision  about  compromise, 
that  must  be  taken  by  people  dealing  with  such  topics. 

In  addition,  the  set  ought  to  overcome  some  of  the  problems 
and  inherent  limitations  of  previously  developed  sets. 

And  of  course,  any  scientific  development  of  a  previous 
theory  or  model  must  yield  better  or  at  least  equivalent 
results . 

A  new  basis  set  has  indeed  been  developed  by  least 
squares  expanding  the  atomic  orbitals  of  atoms,  in  three- 
term  Gaussian  expansions.  The  set  was  optimized  to  give 
the  least  possible  standard  deviations  when  compared  to 
accurate  atomic  orbital  representations.  It  was,  in 
addition,  energetically  optimized  for  best  possible  atomic 
results.  Necessary  adjustments  have  also  been  made 
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before  it  was  applied  to  molecular  systems. 

As  already  mentioned,  the  most  widespread  parallel 
method  used  to  date  is  that  of  ST0-3G.  A  comparison  of 
our  basis  set  with  ST0-3G  on  the  atomic  level  brought 
to  light  the  following  results.  Firstly,  our  basis  set 
invariably  produces  better  atomic  energies.  The  im¬ 
provement  is  significant  as  one  proceeds  to  heavier 
elements.  The  difference  between  our  basis  set  and  ST0-3G 
therefore  increases  as  we  travel  along  the  periods  or  down 
the  groups  of  the  Periodic  Table.  Secondly,  the  atomic 
orbital  energies  produced  by  the  basis  set  developed  here 
significantly  exceed  the  accuracy  of  those  produced  by 
ST0-3G.  The  difference  again  becomes  particularly  worth 
noting  for  the  outer  valence  shells  of  atoms.  The 
importance  of  having  "good"  sets  of  valence  atomic 
orbitals  is  obvious  since  these  shells  are  the  main 
participants  in  molecular  formation.  Still  at  the  atomic 
level,  we  have  demonstrated  the  absolute  inefficiency  of 
STO-SZ  in  treating  transition  metals.  As  expected,  ST0-3G 
should  also  be  judged  inappropriate  for  use  on  third-row 
atoms.  The  superiority  of  the  present  basis  sets  for 
treatment  of  these  atoms  is  unguestionable . 

On  the  molecular  level,  it  has  been  seen  that  the 
results  of  the  present  method  are  once  again  superior 
to  those  of  ST0-3G ,  where  molecular  energies  are  con- 
cerned .  The  binding  energies  and  calculated  eguilibrium 


, 
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geometries  are  at  least  comparable  and  in  many  cases 
better  than  ST0-3G.  Some  exceptions  of  course  do  exist. 
However,  the  resulting  molecular  orbital  energies  are 
consistently  better  than  ST0-3G,  being  closer  to  the  H.F. 
values . 

We  can  thus  safely  conclude  that  the  present  method 
can  successfully  be  applied  to  a  number  of  molecules/ 
giving  very  encouraging  results.  Its  outstanding  advantage 
is  its  efficiency  in  treating  transition  metal  elements. 

It  can  be  argued,  of  course,  that  ST0-3G  is  more  economical 
in  terms  of  computer  time.  However,  the  insignificantly 
small  extra  cost  of  our  method  seems  to  be  more  than 
offset  by  the  quality  of  its  results  and  its  successful 
applicability  to  the  third-period  elements. 
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